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We consider discretized Hamiltonian PDEs associated with a Hamiltonian 
function that can be split into a linear unbounded operator and a regular nonlin- 
ear part. We consider splitting methods associated with this decomposition. 
Using a finite dimensional Birkhoff normal form result, we show the almost 
preservation of the actions of the numerical solution associated with the splitting 
method over arbitrary long time, provided the Sobolev norms of the initial data 
is small enough, and for asymptotically large level of space approximation. This 
result holds under generic non resonance conditions on the frequencies of the lin- 
ear operator and on the step size. We apply this results to nonlinear Schrodinger 
equations as well as the nonlinear wave equation. 



Abstract 



Contents 



1 Introduction 



2 



2 Description of the method 



4 



3 Setting of the problem 

3.1 Hamiltonian formalism . 

3.2 Hypothesis 

3.3 Non resonance condition 



7 

7 
9 
10 



4 Statement of the result and applications 

4.1 Main results 

4.2 Examples 



12 

12 

14 



1 



4.2.1 Schrodinger equation on the torus 

4.2.2 Wave equation on the circle . . . 



14 
16 



5 Proof of the normal form result 



18 



1 Introduction 



In this work, we consider a class of Hamiltonian partial differential equations 
whose Hamiltonian functionals H = Hq + P can be divided into a linear un- 
bounded operator Hq with discrete spectrum and a nonlinear function P having 
a zero of order at least 3 at the origin of the phase space. Typical examples 
are given by the nonlinear wave equation or the nonlinear Schrodinger equa- 
tion on the torus. We consider discretizations of this PDEs and denote by 
= Hq^^ -\- p(^) the corresponding discrete Hamiltonian, where K is a 
discretization parameter. Typically, K denotes a spectral parameter in a collo- 
cation method. 

Amongst all the numerical schemes that can be applied to these Hamiltonian 
PDEs, splitting methods entail many advantages, as they provide symplectic and 
explicit schemes, and can be easily implemented using fast Fourier transform if 
the spectrum of Hq expresses easily in Fourier basis. Generally speaking, a 
splitting scheme is based on the approximation 



associated with the Hamiltonian function Q. For a given time t = nh, n € N, 
the solution starting at some initial value is then approximated by 



The understanding of the long-time behavior of splitting methods for Hamil- 
tonian PDEs is a fundamental ongoing challenge in the field of geometric inte- 
gration, as the classical arguments of backward error analysis (see for instance 
[17]) do not apply in this situation, where the frequencies of the system are ar- 
bitrary large, and where resonances phenomena are known to occur for some 
values of the step size. For example, considering the case of the Schrodinger 
equation on the one dimensional torus, the eigenvalues of -f^Q^'' range from 1 
to and the assumption hK^ << 1 used in the finite dimensional situation 
becomes drastically restrictive in practice. 

Recently, many progresses have been made in the understanding of the long 
time behaviour of numerical methods applied to Hamiltonian PDEs. A first re- 
sult using normal form techniques was given by Dujardin & Faou in [6] for 





(1.2) 
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the case of the hnear Schrodinger equation with small potential. Concerning 
the nonlinear case, results exists by Cohen, Hairer & Lubich, see [8, 9], for 
the wave equation and Gauckler & Lubich, see [10, 11], for the nonlinear 
Schrodinger equation using the technique of modulated Fourier expansion. How- 
ever to be valid these results use non-resonance conditions that are generically 
satisfied only under CFL conditions linking the step-size h and the highest fre- 
quencies of the discretized Hamiltonian PDE. 

Normal form techniques have proven to be one of the most important tools 
for the understanding of the long time behaviour of Hamiltonian PDE (see [1, 
4, 14, 2, 3, 15]). Roughly speaking, the dynamical consequences of such results 
are the following: starting with an small initial value of size e in a Sobolev space 
i?^, then the solution remains small in the same norm over long time, namely 
for time t < Cr£~^ for arbitrary r (with a constant Cr depending on r). Such 
results hold under generic non resonance conditions on the frequencies of the 
underlying linear operator Hq associated with the Hamiltonian PDE, that are 
valid in a wide number of situations (nonlinear Schrodinger equation on a torus 
of dimension d or with Dirichlet boundary conditions, nonlinear wave equation 
with periodic or Dirichlet conditions in dimension 1, Klein Gordon equation on 
spheres or Zoll manifolds.). 

This work is the first of a series of two. 

In this paper, we consider full discretizations of the Hamiltonian PDE, with 
a spectral discretization parameter K that is finite but large. We show that 
under the hypothesis K < £^°' for some constant a depending on the precision 
degree r then the actions of the initial value are almost preserved over a very 
large number of iterations n < Cr£~^ , provided the initial solution is small (of 
order e) in norm. These actions can be interpreted as the oscillatory energies 
corresponding to an eigenvalue of Moreover, the norm of this numerical 

solution remains small for this large number of iterations. 

The method used in this situation is by essence a finite dimensional Birkhoff 
normal form result (explaining why we work here essentially with the norm). 
Using a generic non resonance condition on the step size that turns out to be 
valid for many equations and that is independent on K, we mainly show that 
we can take K asymptotically large without altering the nature of the classical 
result. Our main result is given by Theorem 4.1. 

Roughly speaking, the method consists in applying techniques that are now 
standard in normal form theory, by tracking the dependence in K of the con- 
stants appearing in the estimates. The use of a non resonance condition that is 
independent of K is however crucial, and reflects the infinite dimensional nature 
of the initial continuous problem without space approximation. 

In some sense, the second paper [7] studies the case where K > by 
considering the splitting method where no discretization in space is made (i.e. 
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K = +00). The techniques used involve the abstract framework developed in 
[4, 14, 2]. However, instead of being valid for the (exact) abstract splitting (1.1), 
we have to consider rounded splitting methods of the form 



(1.3) 



where II^^s puts to zero all the frequencies whose weighted energy 
in the Sobolev space H'^ is smaller than a given threshold ry^. Hence, for small r/, 
(1.3) is very close to the exact splitting method (1.1). The good news is that this 
threshold can be taken of the order e"^ , making this projection H^^^ very close to 
the identity, and in any case producing an error that is far beyond the round-off 
error in a computer simulation (particularly for large s). 



Before going on into the precise statements and proofs of this work, we would 
like to give tentative explanations of the restrictions observed in comparison with 
the continuous case. 

The method used in [4] to prove the long-time conservation of Sobolev norms 
and the associated weighted actions for small data is to start from a Hamiltonian 
H = Hq + P depending on an infinite number of variable (^j, r/^), j, G N, and 
for a fixed number r, to construct a Hamiltonian transformation r close to the 
identity, and such that in the new variable, the Hamiltonian can be written 



where Z is a real Hamiltonian depending only on the action /-,• = S^jrjj and R a 
real Hamiltonian having a zero of order r. 

The key for this construction is an induction process with, at each step, the 
resolution of an homological equation of the form 



where G is a given homogeneous polynomial of order n, and where Z, depending 
only on the actions, and x unknown. Assume that G is of the form 

G = Gjk • • • ^jpVki • • • Vkg 

where Gjk is a coefficient, j = {ji, ■ ■ ■ , jp) G and fc = (fci, . . . , kg) S N"?. Then 
it is easy to see that the equation (2.2) can be written 



2 Description of the method 



Ho + Z + R 




{Ho,x] + Z = G 




^U, k)xjk + Zjk = Gjk 



(2.3) 



where 



n{j, k) =ujj^-\ h ujjp -ujki ^j. 
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and where Zj}, and Xjk are unknown coefficients. 

It is clear that for j = k (up to a permutation), we have ^{j, k) = which im- 
poses Zjk = Gjk- When j ^ k (taking into account the permutation), the solu- 
tion of (2.3) relies on a non resonance conditions on the small divisors r2(j, k)~^. 

In [4], Bambusi & Grebert use a non resonance condition of the form 

Vj/fc, \n{j,k)\>jfi{j,kr'^ (2.4) 

where ^(j, k) denotes the third largest integer amongst . . . , They more- 
over show that such a condition is guaranteed in a large number of situations 
(see [4], [14] or [2] for precise results). 

Considering now the splitting method y?^^ o ipp, we see that we cannot work 
directly at the level of the Hamiltonian. To avoid this difficulty, we embed the 
splitting into the family of applications 

[0, 1] 9 A ^ 99^,, o (^^p 

and we derive this expression with respect to A, in order to work in the tangent 
space, where it is much easier to identify real Hamiltonian than unitary flows. 

This explains why we deal here with time-dependent Hamiltonian. Note that 
we do not expand the operator ip^ in powers of h, as this would yields positive 
powers of the unbounded operator Hq appearing in the series. Unless a CFL 
condition is employed, this methods do not give the desired results (and do not 
explain the resonance effects observed for some specific values of h). 

Now, instead of (2.2), the Homological equation appearing for the splitting 
methods is given in a discrete form 

X°^Ho-X + Z = G. (2.5) 

In terms of coefficients, this equations yields 

ie"''"^''''^ -l)Xjk + Zjk = Gjk- 

The main difference with (2.3) is that we have to avoid not only the indices 
{j,k) so that Q(j,k) = 0, but all of those for which h^l{j,k) = 2rmr for some 
(unbounded) integer m. 

In the case of a fully discretized system for which ^ ZjP = for \ j\ > ET, then 
under the CFL-like condition of the form hK^ < C where m depends on the 
growth of the eigenvalues of Hq and C depends on r, then we have \ h^l{j, k)\ < vr, 
and hence 

|gih,t7y,fc) _ ^1 > h^^{^j^-k)-'^ (2.6) 

(2.6) is then a consequence of (2.4). Under this assumption, we can apply the 
same techniques used in [4] and draw the same conclusions. This is the kind of 
assumption made in [8] and [11]. 
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The problem with (2.6) is that it is non generic in h outside the CFL regime. 
For example, in the case of the Schrodinger equation, the frequencies of the 
operator Hq are such that ujj ~ j^. Hence, for large A^, if (ji, . . . , jp, ki, . . . , kg) 
is such that ji = N + 1, ki = N and all the other are of order 1 (N is large here), 
we have il.{j,k) ~ (A^ + 1)^ - A^^ ^ 2N. Hence, 

|gih.ny,fc) _ i| ^ |g2ihjv _ 

cannot be assumed to be greater than /i7/x(j,fc)~" ~ h for all (large) A^. Note 
that a generic hypothesis on h would be here that this small divisor is greater 
than h'~fN~°' for some constants 7 and a. This example shows that we cannot 
control the small divisors ^e'^f^^U^i") _ i| associated with the splitting scheme by 
the third largest integer in the multi index (which is actually of order 1 in this 
case), but by the largest. 

Using a generic condition on h < Hq, we prove in [7] a normal form result and 
show that the flow is conjugated to the flow of a Hamiltonian vector field of the 
form (2.1), but where Z now contains terms depending only on the actions, and 
supplementary terms containing at least two large indices. Here, large means 
greater than e~'^ where a depends on r. 

In the case of a full discretization of the Hamiltonian PDE with a spectral 
discretization parameter K, we thus see that K < £~'^ then the normal form 
term Z actually depends only on the actions, as the high frequencies greater that 
£~'^ are not present. This is essentially the result of this paper. 

In the case where K > £~'^, the normal form result that we obtain can be 
interpreted as follows: the non conservation of the actions can only come from two 
high modes (of order greater than e"'^) interacting together and contaminating 
the whole spectrum. The role of the projection operator H^^^ is to destroy these 
high modes at each step. As we can take rj = e'^, the error induced is very small, 
and in particular, far beyond the round-off error in the numerical simulation. 
This is mainly the result in [7]. 

The differences between the present work and [7] lie in the techniques in- 
volved: In this work, the system considered are large but finite dimensional 
systems, and all the hypothesis made on the nonlinearity can be expressed using 
elementary conditions similar to those used in the finite dimensional case. In [7], 
we study K = +00, which requires much more elaborate technical tools in the 
spirit of [4, 14, 2]. 
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3 Setting of the problem 
3.1 Hamiltonian formalism 

We set N = oi N'^. For a = (ai, . . . , ad) G N, we set 

\a\ = max |aj|. 

i=l,...,d 

Let X G N, and let Mk a finite subset of AA, included in the ball {a G AA| |a| < 
K]. Typically, we can take Mr of the form [-K, . . . , K]'^ C Z'^ or [0, . . . , K]'^ C 
N'^ or a sparse set of the form (see for instance [12, 18]) 

Mk = {a = {ai,...,ad) G 1'^ \ (1 + |ai|) • • • (1 + |arf|) < K} C II^. 

We consider the set of variables {Ca,rib) G C"'^^'^ x C''^^ equipped with the 
symplectic structure 

i dCa/\drja. (3.1) 

We define the set Zk = Mr x {±1}- For j = (a, 5) G Zr, we define |j| = \a\ and 
we denote by j the index (o, —5). 

We then define the variables {zj)j^zK ^ by the formula 



j = (a, 5) G Z;^ 



if <5 = 1, 

7?a if (5 = -1, 



By abuse of notation, we often write 2 = (.^, ry) to denote such an element. 
We set 

II l|2 I |2 

and for any /3 > 0, 

Ba'(p) = {^gC^^|||z|| <p}. 

Note that in the case where K = +00, we set by convention Zr = Z = AAx{ibl} 
and the previous norm defines a Hilbert structure on i'^. We denote by 

iiR-.el^ (c^*-,|| • II ) 

the natural projection. 

Let Ur be a an open set of C^^ . For a function F in C^(Ur,C), we define 
its gradient as 

V^(.) ^ (f ) 



7 



where by definition, we set for j = (a, 6) G Mk x {il}, 



OF 



dF 
dF 

— if 6 = -l. 



Let H(z) be a function defined on Uk- If H is smooth enough, we can associate 
with this function the Hamiltonian vector field Xh{z) defined as 

Xh{z) = JVH{z) 

where J is the symplectic operator induced by the symplectic form (3.1). 
For two functions F and G, the Poisson Bracket is defined as 

^^rx^^ . V- dF dG dF dG 

We say that z G C^^' is real when = for any j G Zk- In this case, 
z = (^, Q for some G C-^-f^'. Further we say that a Hamiltonian function H is 
rea/ if ^/^(z) is real for all real z. 

With a given function H G C°°(Uk, C), we associate the Hamiltonian system 



JVH{z) 



which can be written 



In this situation, we define the flow '^^^{z) associated with the previous system 
(for times t > depending on z G Uk)- Note that if 2: = (^, ^) and H is real, 
the flow {£,^,rf) = ip\j{z), for all time where it is defined, satisfies the relation 
^* = , where ^* is solution of the equation 

dH - 

L = -ij^{(,0, aeATx- (3.3) 
oria 

In this situation, introducing the real variables Pa and qa such that 

^a = ^{Pa + iqa) and = -^{Pa - iQa) , 
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the system (3.3) is equivalent to the system 

Pa = -^{q,p) aeNK 

oqa 

dpa 

where H{q,p) = H{^,^). 

Note that the flow r* = of a real Hamiltonian x defines a symplectic map, 
i.e. satisfies for all time t and all point z where it is defined 

{D,T')U{Dy), = J (3.4) 

where denotes the derivative with respect to the initial conditions. 
The following result is classical: 

Lemma 3.1 Let Uk o,nd Wr- be two domains of C^^' , and let t = ip^ G 
C°°(UK,y^K) be the flow of the real smooth Hamiltonian x- Then for H € 
C°°(>V/<,C), we have 

yz^U XhoAz) = {D,t{z))-^Xh{t{z)). 

Moreover, if H is a real Hamiltonian, H o t is a real Hamiltonian. 

3.2 Hypothesis 

We describe now the hypothesis needed on the Hamiltonian H. 

In the following, we consider an infinite set of frequencies {oJa)a&M satisfying 

VaGAA, \uJa\<C\ar (3.5) 

for some constants C > and m > 0. 

Let lA be an open domain of £^(C^) containing the origin, and let Uk = 
its projection onto C^^. 

We consider the collection of Hamiltonian functions 

ijW = nf^ + pW, K > 0, (3.6) 

with 

Hf^ = uJaiz) 

a&NK 

where for all a G Mr-, 

Ia{z) = ^aVa (3.7) 
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are the actions associated with a £ Mk- Note that if z = (C, then Ia{z) = iCaP- 
We moreover assume that the fmictions P(^) G C°°{Uk,€.) are real, of order 
at least 3, and satisfy the following: For all ^ > 1, there exists constants C{i) > 
and > such that for all K > 1, (ji, • • • , j^) E Zj^ and z G ^i^-, the following 
estimate holds: 

< C{e)K'^^^\ (3.8) 
The Hamiltonian system (3.2) can hence be written 

I Va 

Denoting by (/5q the exact flow of a Hamiltonian flow, splitting methods are based 
on the approximation 

^0 

for a small time step h > 0. Note that in this case, the exact flow of Hq^^ is 
explicit and given by 

while the calculation of y^^(K) requires the solution of an ordinary differential 
equation, whose solution is often given explicitely (see the examples below). 

The goal of this paper is the study of the long-time behavior of the numerical 
solution given by (1.2) for large number n of iterations. 

Remark 3.2 Note that no hypothesis is made here concerning the preservation 
of the norm by the flow of (3.9). 

3.3 Non resonance condition 

In the following, for j = {ji, ■ ■ ■ , jr) G Z^ with r > 1, we use the notation 

~ ^jl ' ' ' ^jr ■ 

Moreover, for j = {ji,...,jr) G Z'k with ji = {ai,5i) £ Mk x {±1} for i = 
1, . . . , r, we set 

j = (ji,---Jr) with J. = (ai,-(5i), i = l,...,r. 



dzj, ■ ■ ■ dzj^ 



orja 

Oc,a 



(3.9) 
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and we define 

= SiUJai H h drUJar- 

We say that j G Z'^ depends only of the action and we write j G if r is even 
and if we can write (up to a permutation of the indexes) 

Vz = l,...r/2, ji = {ai,l), and ij+r/2 = (a*, -1) 
for some Oj G Nk- Note that in this situation, 

^3 ~ ^ji ' ' ' ^jr ~ CaiVai ■ ■ ■ ■Ca^/a^a^/a 

where for all a £ Mk, Ia{z) denote the actions associated with a (see (3.7)). For 
odd r, is the empty set. 

We will assume now that the step size h satisfies the following property: 

Hypothesis 3.3 For all r G N, there exist constants 7* and a* such that for all 
K G N*, 

(ii,...,jV,)GZ]:A-45, =^ |l_e^'^^(i)|>-|rl. (3.10) 

The following Lemma 3.5 shows that condition (3.10) is generic in the sense 
that it is satisfied for a large set oi h < Hq (and in particular independently of 
K), provided that the frequencies uJa satisfy a non resonance condition that we 
state now (see [17, 19] for similar statements): 

Hypothesis 3.4 For all r G N, there exist constants j{r) and a{r) such that 
G N*, 

in, . . . ,jr) e z^j,\a'k =^ > ^. (3.11) 

In the next section, we will check that condition in different concrete cases. 

Lemma 3.5 Assume that Hypothesis 3.4 holds, and let Hq and r be given num- 
bers. Let 7 and a be such that (3.11) holds and assume that 7* < (2/7r)7, 
a* > a + ma + r with a > 1 and m the constant appearing in (3.5), then we have 

* 

meas{ /i < /iq | /i does not satisfy (3.10) } < C — 

7 

where C depends on a and r. As a consequence the set 

Z^ho) = {h < hQ\h satisfies Hypothesis 3.3} 
is a dense open subset of (0, /iq). 

The proof of this lemma is given in [7, Lemma 4.6]. 
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4 Statement of the result and applications 
4.1 Main results 

Theorem 4.1 Assume that p(^) and h < Hq satisfy the previous hypothesis. 
Let r G N* be fixed. There exist constants a, C and Eq depending only on r, /iq 
and the constants (3{l) and C{£), i = 0,...r in (3.8), such that the following 
holds: For all e < Eq and K < e"'^ , and for all real such that 

if we define 

z"=(^^(,,o^^(,,)"(z°) (4.1) 
then for all n, is still real, and moreover 

\\z^ < 2e for n < (4.2) 

and 

yaGMK, \Ia{zn - Ia{z^)\ < Ce^^^ for n<^ (4.3) 

The proof of this resuh rehes on the following Birkhoff normal form result, 
whose proof is postponed to Section 5: 

Theorem 4.2 Assume that that p(^) and h < ho satisfy hypothesis (3.8) and 
(3.10). Let r G N* be fixed. Then there exists constants [i and C depending on r, 
ho, and C{i), i = 0, . . . r in (3.8) and a canonical transformation tk from 
Bk{p) iiT'to BkC^p) with p = {CK)~^ satisfying for all z G Bk{p), 

\\tk{z) - z\\ < {CK)f^\\z\\^ and \\t^\z) - z\\ < {CKf\\zf , (4.4) 

satisfying the following result: For all z £ Bk{p), 

'^K^ ° "ftik) ° fp(K) o rxiz) = (p^K) o iPk{z) 

where ipK satisfies: 

• ipxiz) is real if z is real, 

• For all z G Bk{p), 

Uk{z)-z\\ <{CKf\\zf . (4.5) 

• For all z G Bk{p), 

\Ia{ll^K{z)) - La{z)\ < {CKf\\z\C' . (4.6) 
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Proof of Theorem 4.1. First, let us note that as the Hamiltonian functions 
Hf^ and are real Hamiltonians, it is clear that there exist £ C-'^ such 

that for all n, we have = (^",^'^), that is is real. 
Let (3 given by Theorem 4.2 and let a = 1/(2/3). We have for K < s~'^, 

{CKf < 

Let tk be defined by Theorem 4.2, and let = r^^(z"). Using the property of 
tk, we see that is real, i.e. we have = for all n. By definition, we 

have 

Vn>0, y"+^ = (c^^K)oV'x)(y"). (4.7) 

Using the fact that K < e~°" and (4.4), the transformation tk in the previous 
Theorem satisfies the following: For all z such that ||z|| < 2e, 

\\t^\z) - z\\ < C^e-i/2||zf < 4C^e3/2 < 1^ (4^8) 

provided Eq is sufficiently small. Hence we have \\y^\\ = \\t^^{z'^)\\ < |e. 

Note that we have p = {CK)~^ > C~^e^^'^ > 2e provided that eq is small 
enough. Using (4.5) we get that as long as < 2e, we have 

< ||y"|| +{CKf\\y''f < ||y"|| + 2^CV-^/2 
By induction, we thus see that for 

we have \\y"'\\ < jE < 2e. Assuming that Eq is such that 2~^~^C~^eI^'^ < 1, 
this shows that for n < e^~'' we have lly*^!! < |e. Using (4.4) and an inequality 
similar to (4.8), we conclude that 

||z"|| < 2e, for n < 

II II - > - ^r-l 

which yields to (4.2). 

Now using (4.6) and the fact that ||y"|| < 2e we see that for n < e^~'^ we have 
yaeAfK, |/a(y"+') - Ia{yn\ < 2^+1^7^+1/2 

whence 

Va G A/-K, \Ia{yn - Ia{y^)\ < 2^+iC^ne^+i/2 
Now we have for all a E Mk 

\uyn - ia{zn\ = \\a' - \a'\ = \\a - x \\a + \ai 

whence 

\iaiyn-Un\ < \C-C\i\\y''\\+\\z^) < WTKiyn-y'^W (lly'^ll + 11^^11 )• 



13 



Using (4.4) we see that for all n < ^' and all a G Mk, 

||Ti^(i/") <4C^e3/2. 

and hence, as \\z^\\ < 2e, 

\Ia{yn-Ia{zn\<8Cf'e'/^. 

Using (4.8), we thus see that 

Vn < e^~'- Va G Mk, |/a(z") - Ia{z^)\ < 2''+^C^{£^/^ + ne'+^Z^) 
and this easily gives the result. ■ 

4.2 Examples 

In this section we present two examples, other examples like the Klein Gordon 
equation on the sphere (in the spirit of [3] ) or the nonlinear Schrodinger operator 
with harmonic potential (in the spirit of [15]) could also be considered with these 
technics. 

4.2.1 Schrodinger equation on the torus 

We first consider nonlinear Schrodinger equations of the form 

idtip = -A^ + Vi^Tp + d2g{^,ti^), xeT"^ (4.9) 

where V G C°°{T'^, M), g G C°°{l/(, C) where ^ is a neighborhood of the origin in 
C^. We assume that g{u,u) G M, and that g{u,u) = 0(1^^). The corresponding 
Hamiltonian functional is given by 

H{'iP,i,)= [ \V'ilj\'^ + ij{V + g{i;,i;)dx. 

Let 4>a{x) = e*°'^, a G Z"^ be the Fourier basis on ^^(T'^). With the notation 
1 ^ _ _ ^ ^ / 1 \d/2 



/ 1 Nd/2 _ / 1 Nd/2 

^ = \2t^) ia4'a{x) and V = [^) 2^ ??a0a(2;) 



the (abstract) Hamiltonian associated with the equation (4.9) can be formally 
written 

H{i, r,)=J2 ^aiaria + P{i: T,). (4.10) 

Here = |«P + are the eigenvalues of the operator 

tp I— > —Alp + V -kip, 
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and we see that uja satisfy (3.5) with m = 2. Moreover, the nonhnearity function 
P(^, 77) posesses a zero of order 3 at the origin. In this situation, it can be shown 
that the Hypothesis 3.3 is fulfihed for a large set of potential V (see [4] or [14]). 

Following [10], a space discretization of this equation using spectral colloca- 
tion methods yields a problem of the form (3.6) with 

Nk = [-K,...,K-lY 

and, with 

/ 1 Nd/2 , I xd/2 _ 

UK = [i^j 2^ ia(t>a{x) and VK = y—j 2^ r]a(l}a{x) (4.11) 
the nonlinearity reads 

P^^'Ki^r^) = I Q[9{uK,VK))dx 
where, for a function ij; = {^Y^'^ Y.a&d i'a(t)a{x) 

2W= V'a+2i^fe)</'a(x) 

is the collocation operator associated with the points Xa = ^ € T*^, a G J^k- It 
is easy to verify that p(^) satisfies (3.8) for some constants C{£) depending on 
g and depending on g and the dimension d. 

Note that starting from a real initial value u^{x) (see (4.11)) this system 
reduces to solving the system of ordinary differential equation 



\/a^Na i-^UK{Xa,t) = J^2K^J^2K'^i< i^a, t) + d2g{uK{Xa, t) , UK {Xa,t)) 

where is the matrix (u>a)a£j\fK ^''^^ ^2K the Fourier transform associated 
with Nk- Note that in this case, the numerical solution (1.2) is easily im- 
plemented: The linear part is diagonal and can be solved explicitely in the 
Fourier space, while the non-linear part is an ordinary differential equation 
with fixed parameter Xa at each step. If moreover g{u,u) = G(|up) for some 
real function G then the solution of the nonlinear part is given explicitely by 
V^{K){u) = ex.p{—2ihG' {\u\'^)u using the fact that |np is constant for a fixed 
point Xa- 

For high dimension d, the previous discretization is usually replaced by a 
discretization on sparse grid, i.e. with 

AAx = { a = (ai, . . . , ad) G I (1 + |ai|) • • • (1 + \ad\) <K}c Z"'. 
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As explained in [18, Chap III.l], methods exist to write the corresponding system 
under the symplectic form (3.6), upon a possible loss in the approximation prop- 
erties of the exact solution of (4.9) by the solution of the discretized Hamiltonian 
H^^\ Note that this does not influence the long time results proven here: In 
some sense we do not impose the nonlinearity P^^\z) to approximate an exact 
nonlinearity P{z). 

We give first a numerical illustration of resonance effects. We consider the 
equation 

in the one dimensional torus T"*^, with initial value 

o 



— cos(x) 

Note that this problem is equivalent to solving (4.9) with a small initial value 
of order e. We take e = 0.1, V with Fourier coefficients Va = 2/(10 + 2a^) and 
K = 200 (i.e. 400 collocation points). In Figure 1, we plot the actions of the 
numerical solution given by the Lie splitting algorithm (1.2) in logarithmic scale. 
In the right we use the resonant stepsize h = 27:/ {uj — wi) ~ 0.17459 .... In the 
left we plot the same result but with the non resonant stepsize h = 0.174. 




Figure 1: Plot of the actions for non-resonant and resonant step size. 



In Figure 2, we show the long time almost conservation of the action in the 
case where /i = 0.1 (non resonant), and e = 0.1 and e = 0.01 after 10^ iterations. 



4.2.2 Wave equation on the circle 

We consider the wave equation on the circle 
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9 
s 

I "° 

9 



Figure 2: Conservation of the actions for e = 0.1 (left) and e = 0.01 (right). 



where m is a non negative real constant and g a smooth real valued function. 
Introducing the variable v = ut, the corresponding Hamiltonian can be written 

H{u,v) = I -{v'^ + ul + mu^) + G{u)dx, 
J J 2 

where G is such that duG = g. Let A := {—dxx + m)^/^, and define the variables 
{P,q) by 

q ■= A^^'^u, and p = A'^l'^v. 
Then the Hamiltonian can be written 

i/=i((Ap,p)^2 + (Ag,g)^2)+ / G(A-V2g)dx. 

Let tjja = \/\a\'^ + m, a G N =: be the eigenvalues of the operator A, and 
the associated eigenfunctions. Plugging the decompositions 

= '^<la(t}a{x) and p{x) = '^Pa(l)a{x) 

aSN aGN 

into the Hamiltonian functional, we see that it takes the form 

aGN 

where P is a function of the variables pa and qa- Using the complex coordinates 

1 1 
ia = ^{qa + iPa) and rja = -j={qa - iPa) 
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the Hamiltonian function can be written under the form (4.10) with a nonlinearity 
depending on G. As in the previous case, it can be shown that the condition 
(3.11) is fulfilled for a set of constant m of full measure (see [4, 2]). A collocation 
discretization on equidistant points of [0, 27r] yields the same discretization as 
previously (with d = 1). 

In this situation, the symmetric Strang splitting scheme 

h/2 h h/2 

corresponds to the Deuflhard's method [16]. If moreover we consider the Hamil- 
tonian 

h'^^\z) = Hf\z) + P^^\^{h^l)z) 

where VL is the matrix with elements uja, cl ^ -^K, and <I>(x) a smooth function that 
is real, bounded and such that $(0) = 1, then the splitting schemes associated 
with this decomposition coincide with the symplectic mollified impulse methods 
(see [17, Chap. XIII] and [9]). 

5 Proof of the normal form result 

The rest of the paper consists in proving Theorem 4.2. 

In the following, we denote by %• the set of polynomial of order r on C^^' (for 
sake of simplicity, we do note write the dependance in K in the notation 7^). If 



is an element of we set 



\Q\^ = max max \Q-j\. 



r £=0,...,rje2^ 

If moreover Q £ C{[0, !],%■) we set 

IIQII^ = max |Q(A)|^. 

Ag[0,1] 

Using the assumptions on P^^^ , we can write a Taylor expansion of P around 

0, 

r 

P(^){z) =P, + Q^ = Y^ PjZj + Qr{z) 

where 

\Pj\< CKf^° 
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where C and (3q depend on I3{1) and C(£), £ = 0, . . . , r in (3.8). 

Notice that Qr{z) G C~(C^^,C) admits of zero of order r + 1 and satisfies 

||Xq,(z)|| <CK/^"||z|r 

for z G Uk, provided (3q = (5o{r,d) is large enough. 

Before giving the proof of Theorem 4.1, we give easy results on the flow of 
non autonomous polynomials Hamiltonian. 

Lemma 5.1 Let k>l and let P{\) G C([0, 1], 7fc_|_i) he a homogeneous polyno- 
mial of order k + \ depending on X G [0, 1]. Then 

(i) There exists a constant C depending on k such that for all z G and all 

A G [0, 1], we have 



\P{X,z)\ < 



I iik+i 
\z\\ 



fe+i 



(a) There exists a constant C depending on k such that for any z G and 
all A G [0,1], 



\z\ 



,k 



fc+i 



Moreover, Let ki and /c2 two fixed integers. Let P and Q two homogeneous 
polynomials of degree ki + I and /c2 + 1 such that P G C([0, l],7fc^-|_i) and Q G 
C([0, l],7fc2+i). Then {P,Q} G C([0, 1], T^j+fc^) and we have 

\\{P,Q}\\^^ <C||P|| IIQII 

for some constant C depending on ki and k2- 
Proof. We have 

where we have set for j = {ji, . . . G 2^^, 

Using \zj\ < ll^ll we easily obtain (i) using ^Zk < + 1)"^. The second 
statement is proven similarly. The estimate on the Poisson brackets is trivial. ■ 



Lemma 5.2 Let r > 3, 

r 

Q(A,z) = ^ E QiW^i 
^=3 je^l. 
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be an element of €{[0,1],%). Let ^q(^x) A^'^ associated with the nan 

autonomous real Hamiltonian (5(A). Then there exist a constant Cr depending 
on r such that 

p<mi{l/2,CrK-^'-\\Q\\-') =^ VAg [0,1], C i?x(2p). 

(5.1) 

Moreover, if F{X) £ C{[0,1],C°° {Bk{p),C)) has a zero of order r at the origin, 
then F{\) o ^q(^x) ^ '^f order r at the origin in Bk{p)- 

Proof. Let z'^ = Lp^^y^-^{z^) . Using the estimates of the previous lemma, we 
have 

A||^A||2 ^ 2(z\Xq(,)(z^)) 

for some constant Cr depending on r. Hence, as long as H^'^H < 1, we have 

By a standard comparison argument, we easily get that for G Bk{p) we have 

VA G [0,1], llz^ll < 2||z°|| . 
This shows (5.1) and the rest follows. ■ 

We give now the general strategy of the proof of the normal form Theorem 4.2, 
showing in particular the need of working with non autonomous Hamiltonians 
and of considering the non resonance condition (3.10). 

We consider a fixed step size h satisfying (3.10). As in this section K will be 

considered as fixed, we denote shortly by P and hI^^ by Hq. We consider 
the propagator 

h h h I 

We embed this application into the family of applications 

v'ho^^hP, AG [0,1]. 

Formally, we would like to find a real Hamiltonian x = x(A) and a real Hamilto- 
nian under normal form Z = Z{X) and such that 

V A G [0, 1] (/.^^ o v9^p o (^J(^) = (^J(^) o o . (5.2) 

Let G C^^' and z'*' = ip^^ o ip^p o cp^^^^^z^). Deriving the previous equation 
with respect to A yields 
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Using Lemma 3.1 that remains obviously valid for non autonomous Hamiltonian, 
we thus have 

= Xa{x){z ) 

where ^(A) it the time dependent real Hamiltonian given by 

A{x) = hP o (^-; + x{\) o ipi^ o 

Using the same calculations for the right-hand side, (5.2) is formally equivalent 
to the following equation (up to an integration constant) 

V A G [0, 1] hP o c^-J + x(A) o o c^-J = x{X) + hZ{\) o ^-^^^ o c^^J. (5.3) 

which is equivalent to 

VAg[0,1] xWoA,-xWo^lp = hP-hZ{\)o^-^^^. (5.4) 

In the following, we will solve this equation in x(A) and Z{\) with a remainder 
term of order r + 1 in z. So instead of (5.4), we will solve the equation 

V A G [0, 1] x(A) o (^^^ - x(A) o (^-^ = hP- {hZ{\) + R{\)) o ^-J^^y (5.5) 

where the unknown are x(A), and Z{\) are polynomials of order r, with Z under 
normal form, and where -R(A) possesses a zero of order r + 1 at the origin. 
In the following, we formally write 

r r 

x(A) = ^X[e\W = XjWzj 

e=3 £=3 j^z\^ 

and 

r r 

where here the coefficients Zj{X) are unknown and where we denote by X[^](A) 
and Z[gj{X) the homogeneous part of degree £ in the polynomials x(A) and Z{X). 
Identifying the coefficients of degree £ < r in equation (5.5), we obtain 

Xm(A) o ifH^^ - X[e]W = hP[e] - ^%](A) + /iG[^] (A; x*, -P*, no- 
where G is a real Hamiltonian homogeneous of degree i depending on the poly- 
nomials X[k]: P[k] and Zffc] for k < L In particular, its coefficients are polynomial 
of order < ^ of the coefficients Xj; Pj and Zj for j G Z^, k < i. 
Writing down the coefficients, this equation is equivalent to 

yj eZl^ (e'^^ij) _ 1)^^. = hPj - hZj + hGj 

and hence we see that the key is to control the small divisors e^^^^^^ — 1 to solve 
these equations recursively. 
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Lemma 5.3 Letxi^) be an element of €{[0,1],%.). Letr^X) := ip^^^^ be the flow 
associated with the non autonomous real Hamiltonian x{^)- Let g G C([0, 1],7^), 
then we can write for all do G [0, 1], 

5(0-0) o T(cro) = giao) 

+ Z] / ■■■ / (Ad^K)0---oAd^(<^,)C/(cro))do-i---dcrfc + ii(cro) (5.6) 
where by definition Adp(Q) = {Q,P} 

-R(o"o) = j ■■■J (Ad^(<^^) o • • • o Ad^(^^)5r((To)) o r(cr,.) ddi • • • dcr^. (5.7) 

Each term in the sum in Eqn. (5.6) belongs (at least) to the space C{[0, 1],%^.). 
The term R{ao) defines an element of C{[0,1],C°° {C^'^ ,C)) and has a zero of 
order at least r + 1 at the origin. 

The proof of this lemma is given in [7] . 

As mentioned previously, for a given polynomial x £ C([0, 1], 7^) with r > 3, 
we use the following notation 

r r 

where X[£](A) G C([0, l],?^) is a homogeneous polynomial of degree I. 

Proposition 5.4 Lei x(A) be an element 0/ C([0, 1], 7^). Let ^\{X) /flow 
associated with the non autonomous real Hamiltonian x(A). Let g € C([0, 1],7^), 
t/ien we can write for all A G [0, 1], 

5(A)o^^(,)=5W(A)+rM(A) 

w/iere 

• ^^(A) eC{[0, !],%.). Moreover, if we write 

r 

e=3 

where S^i^ (A) is a homogeneous polynomial of degree I, then we have for all 
£ = 3,...,r, 
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where G[^](A; x*)5*) is a homogeneous polynomial depending on A and the 
coefficients Sj are polynomials of order < I of the coefficients appearing in 
the decomposition of g and x- Moreover, we have 

\\Gi4^;x.,g*)\\ < (i + EII5[hII')(^+EII^mI0- ^^'^^ 

m=3 m=3 

• T('")(A) G C{[0,l\,C'^{C^^ ,C)) has a zero of order at least r + 1 at the 
origin and satisfies for all z £ Bk{^/'^), 

where 

r r 

Crix*,9.) <c{l+Y. ll5MllrJ(l+ E II^mIItJ (5-10) 

m=3 m=3 

with C depending on r. 

Proof. Using the previous lemma, we define S^'''^ as the polynomial part of 
degree less or equal to r in the expression (5.6): this polynomial part may be com- 
puted iteratively, from the homogeneity degree 3 to r. Actually, every Poisson 
bracket appearing in (5.6) is taken with a polynomial xi^k)-, which decomposes 
into homogeneous polynomials with degree 3 at least. The terms appearing in 
the sum in (5.6) hence have an increasing valuation, and this allows the iterative 
computation. The remainder terms, together with the term R{X) in (5.7), define 
the term T^*^') (which is an element of C([0, 1], 72^)). The properties of S'^^''(A) 
and T^"^\\) are then easily shown using Lemma 5.1. ■ 

The next result (Proposition 5.6 below) yields the construction of the normal 
form term ijjx of Theorem 4.2. 

Definition 5.5 A polynomial Z on is said to he in normal form if we can 
write it 

r 

where Af^ is defined in the beginning of Subsection 3.3. 

Proposition 5.6 Assume that H := H^^^ satisfies (3.6) with P := ful- 
filling (3.8) and assume that h < Hq satisfies the hypothesis (3.10). Then there 
exist 
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a polynomial x G C([0, 1], 7^^ 



e=3 e=3 j(zzi^ 

• a polynomial Z £ C{[0, !],%■) 

r r 

e=3 e=3 j^A'j, 

in normal form, 

• a function R{\) € C{[0,1],C°° {Bk{p),C)) with p < coK~^ for some con- 
stant Co > and (3 > 1 depending on r and d, and having a zero of order 
at least r + \ at the origin 

such that the following equation holds: 

VAG[0,1] xWoA.-x{\)o^h^ = hP-{hZ{\)+R{\))o^-^^^. (5.11) 
Furthermore there exists a constant Cq depending on r and d such that 

11x11^^ + \\Z\\^^ < CoK^ 
and such that for all p < coK~^ and all z G Bk{p), we have 
VAg[0,1], \\XR(^^-^{z)\\<CoKf\\: 



I T 
\z\\ 



Proof. Identifying the coefficients of degree £ < r in the equation (5.11), we 
get 

X[e\ ° ^Ho - X[e] = hP[i] - hZ[e] + hGii]{x*, P*,Z^). 

where G is a real Hamiltonian homogeneous of degree i depending on the poly- 
nomials X[k]j P[k] Slid for k < I. In particular, its coefficients are polynomial 
of order < ^ of the coefficients Xj; Pj ^^'^ for j G Z'^, k < I and satisfy 
bounds like (5.9). Writing down the coefficients, this equation is equivalent to 

yj^Ir (e*^^(^') - l)xj = hPj - hZj + hGj. 

We solve this equation by setting 

Zj = Pj + Gj and Xj = for j G A^^ 

and 

Zj = and Xj = -j^^^^jy—^iPj + Gj) for j Mx- 
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Using (3.10) and the result of Proposition 5.4 we get the claimed bound for some 
P depending on r. 

To define R, we simply define it by the equation (5.5). By construction and the 
assumption on P = P(^), and using bounds of the form (5.10), it is easy to show 
that it satisfies the hypothesis. ■ 

Proof of Theorem 4.2. Integrating the equation (5.5) in A, it is clear that 
the following equation holds: 

V A G [0, 1] o ip^p o ip^^^^ = (^J(^) o ip%^ o ip^z{\)+R(X) ■ 

Note that using Proposition 5.4 and (5.1) we show that for z G Bk{p) with 
p = cK~^ we have 

||9'J(A)(^)-^II <CK%f. 
This implies in particular that 

\\4 < y'^^x)i^)\\+CK~4z\\ 

For K sufficiently large, this shows that ^^(^x) invertible and send Bxip) to 
Bk{2p). Moreover, we have the estimate, for all A G [0, 1], 

\\{p'^^,^)-\z)-z\\<CK%f. 

We then define tk = V^^(a) ^^^^ ^pK = VhZ(\)+R{x) ^^'^ verify that these applica- 
tions satisfy the condition of the theorem for suitable constant C and /3. ■ 
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